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Abstract 

The anomalous part of the NMR relaxation rate of copper nuclei in the normal state of copper- 
oxide metals is calculated using the orbital magnetic parts of the fluctuations derived in a recent 
theory to explain the long wavelength transport anomalies. Oxygen and Yttrium reside on lattice 
sites at which the anomalous contribution is absent at all hole densities. The frequency, momentum 
dependence, and the form-factor of the fluctuations is predicted, which is verifiable by inelastic 
neutron scattering experiments. 
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Of all the various anomalies in the normal state of copper-oxide metals, the most extraordinary 
are the nuclear magnetic relaxation 2 rates T^ 1 . These rates probe local fluctuations and are a very 
detailed and stringent test of microscopic theories. The anomalies are essentially identical in the 
best studied materials, YBa2Cu306+ x with x « 0.9 (12306.9) and La2-<5Sr,sCu04 with 5 « 0.15 i.e. 
near the composition for the highest T c . The local magnetic fluctuations in a metal are usually of 
the form: Imx(w) ~ leading to T^ 1 ~ T. The oxygen T 1 1 in the Cu-0 metals are consistent 

with this behavior, 3 as is that of Y nuclei in (123). 4 The Cu-nuclear relaxation rate, on the other 
hand, is 5 to a first approximation, CuT^" 1 ~ a + bT with a ~ 6T C . The (nearly) constant part 
in Cu Ti 1 sug gests that the local magnetic fluctuations in Cu have a scale-invariant form at low 
energies: x"(<jj) ~ ^/T, 6 unlike Fermi-liquids and as if the metals lie near a T = critical point. 
This is remarkable enough. But it is even more remarkable that the local electronic fluctuations on 
the oxygen do not share the fluctuations on Cu even though the electronic wavefunctions on the 
two ions are well hybridized. 

One idea discussed extensively 7 to explain the observations is that the materials are near an 
antiferromagnetic instability of the Cu magnetic moments with a magnetic correlation length £ _1 ~ 
T. Since the oxygen ions sit half way between the coppers, they do not see the antiferromagnetic 
fluctuations. Direct measurements 8 of x"(ci> w ) by inelastic neutron scattering show this idea to 
be untenable. In 12306.9 the magnetic correlation length is 8 less than a lattice constant at low 
energies and temperature independent. In Lai.ssSr 15CUO4, magnetic correlations at low energies 
are indeed observed. 9 But the correlations are at an incommensurate wave-vector so that the 
fluctuations do not cancel on the oxygen ions. Calculations 10 of T^f 1 using the measured x"(q,u>) 
give a temperature dependence which is quite different from the experiments for both °Tjf 1 and 
Cu^p-i Moreover, (°TiTK) _1 « (1.4±0.2) sec -1 deg -1 for all measured temperatures in the normal 
phase for 12306+ x for several different x and for Lai.ssSr. 15CUO4. 11 Here K is the measured oxygen 
Knight shift. This is in agreement with the value obtained 12 from the conventional Korringa Law, in 
terms of the deduced hyperfine coupling constants. The temperature independence of this quantity, 
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as well as its near invariability from one compound to another is inexplicable by AFM correlations 
whose q- vector and magnitude change from compound to compound and with x. Some more robust 
symmetry is called for so that the fluctuations seem by Cu nuclei are absent at the oxygen sites. 
Clearly, something more subtle is at work. 

Recently a model for copper-oxide metals 13 ' 14 has been systematically investigated. 14 A T = 
critical point has been discovered at a unique composition x = x c , at which the fluctuations lead 
to marginal Fermi-liquid self energy: ImE(w) ~ max(u;,T) sgnu;, and long wavelength transport 
properties in agreement with experiment. Here these fluctuations are applied to calculate features 
of TY 1 . 

For reasons discussed elsewhere, 13,14 the model includes dynamical degrees of freedom both on 
Cu and on O ions and includes long-range interactions besides a strong short-range repulsion on 
Cu ions to exclude the Cu 3+ state. For long-range interaction energies comparable to or larger 
than the Cu-0 bonding - antibonding splitting, the model has a critical point at T = 0, and hole 
density x = x c . Near this point (intracell Cu-O) current fluctuations have a propagator 14 

D(q, u) = (1/W B )( W7 + In^c/T)}- 1 + Aq 2 + G(x, T))" 1 , (1) 

Here W s is the spectral weight of the fluctuations and j(u, T, G), the damping of this unconserved 
mode is rs max(|u;|, T, G, t^ 1 ). Tq 1 is the order of the elastic rate for single particle scattering 
from impurities, obtainable, from resistivity measurements. The "mass" of the collective mode, 
G(x, T) = go|x — x c | + T/To, where go is 0(1) and To is the order of the bandwidth. 

We now consider the orbital magnetic contribution of these current modes to the magnetic 
fluctuations: x(q, uj). First, consider the local fluctuations around any particular ion with the Cu- 
O metal treated in the one-electron approximation. The local hyperfme Hamiltonian for a nuclear 
spin I chosen to be at the origin is 

H hf = I • M = I ■ M kq (0)c+ +q/2 c k _ q/2 , (2) 

kq 
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where M kq (0) is the matrix element of the magnetization at the origin between Bloch waves: 15 

|k >= e ik r n k (r) (3) 



The contributions to M relevant to us are 

M orb = i^ = i^^^- (4) 

Consider the matrix element 

(k+||M orb |k-q/2) = ^e-^X+q/aW + u k _ q/2 (r)). (5) 

The Bloch- functions u k (r) are in general very complicated. We will content ourselves with 
general properties of (5) remembering that the relaxation rate sums over all k and q. To estimate 
(5) consider u k at some symmetry points where we know its properties, for example the non-bonding 
point k = (0, 0). At this point uo(r) can be written as 

-^£ai, Q (0)<Kr-R iiQ ), (6) 

where i sums over the cells and a the atoms in the cell: </>(r — R\, a ) has the symmetry of the d x 2_ y 2 
orbitals at the Cu-site and has the symmetry of the p x and p y orbitals, appropriately, at the two 
oxygen sites per unit cell. The phase factors ai jQ (0) are ±1, phased to produce the non-bonding 
configuration at k = 0. We can verify by "k.p perturbation" to determine the phase factors ai jQ (k) 
for k about this point that the qualitative results remain the same at arbitrary points in the zone. 
So we can use (6) to find the qualitative properties of (5). 

Inserting (6) in (5), we use the rapid fall-off of <j)(r — Ri Q ) with |r — Ri Q | to approximate (6) by 
a sum of two terms: (a) (i, a) at the origin (i.e. the position of the nucleus whose relaxation rate is 
being evaluated) on both u(r) factors in (5); (b) (i, a) at the origin in one u(r) and (i, a) = nearest 
neighbors of the atom at the origin and vice-versa. In higher order terms the sums over (i, a) may 
be replaced by an integral over space. This is then just the fluctuating electromagnetic field at 
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the origin due to long wavelength orbital currents. Such effects have been evaluated and formed 
to contribute O(10 -4 ) to the relaxation rates compared to the on-site relaxation terms (say, from 
contact hyperfine interactions) and may be ignored. 

The local orbital fluctuation term (a) in (6) is zero at both Cu and oxygen if as in Cu-0 metals 
non degenerate orbitals reside at these sites. Consider now the next term (b). The second 17 term 
in (5) yields a contribution linear in q given by 

(k + §|M orb |k - S) ~ fc^iq . ( r< P(r - R n )^(r)), (7) 
z z n r 

where n sums over the nearest neighbors. Take a Cu-site at the origin; then cf>(i) has d x 2_ y 2 
symmetry so that L z <j)(r) ~ d xy . cf)(r — R n ) have p x or p y symmetries so that XInl " r( ?H r ~~ Rn) 
for a general q also has a component with d xy symmetry. Finite matrix element of L z therefore 
exist. The construction of a d xy symmetry around a Cu ion by appropriately phasing the oxygen 
p x and p y orbitals is illustrated in Fig. 1. Semi-classically, orbital relaxation requires a current 
fluctuation circulating around the nuclei. This is provided by the loop formed by the four oxygens 
around a given Cu. These produce magnetization fluctuations in the z-direction which contribute 
to the nuclear relaxation rate TJ~2 with field applied in the Cu-0 planes. 

The same mechanism also contributes to Cu T^| 1 with field applied normal to the Cu-0 plane. 
Since L Xjy d x 2„ y 2 ~ d xz ,d yz this occurs through the relative phasing of the oxygen p Xjy orbitals in 
the plane and the p z orbital of the apical oxygens with a q in the plane and in general with the 
apical oxygens mutually out of phase. 

Since L z p x ~ p y , etc, orbital relaxation of oxygen nuclei requires a fluctuation of p x>y form to be 
created by phasing the neighbors. For the magnetic fluctuation at the oxygen nuclei, consider (6) 
with oxygen sites at the origin. Oxygen is linearly coordinated by copper ions, and no contribution 
of O(q) exist from the nearest neighbor terms, because J2n Q ' r( P( T ~ R-n) cannot be phased to yield 
a p-symmetry at the oxygen at the origin. This is true to any order in q. Semiclassically, no simple 
closed loop current fluctuation can be created around oxygen unlike the case of copper illustrated 
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in Fig. (1). Similarly, Y is an s-state ion sitting between Cu-0 planes and no contributions of the 
type we are considering occur for it. 

Magnetic fluctuations, which because of the symmetry of the lattice, contributing only to re- 
laxation of Cu-nuclear spins therefore do exist. We now consider their frequency and temperature 
dependence. We use the fact that the operators in (4) project on to operators on Cu and O-orbitals 
which participate in the current fluctuations of Eq. (3). The magnetization fluctuations then couple 
to the current fluctuations through a factor proportional to q, as in Eq. (7). We therefore have an 
anomalous contribution to the relaxation at Cu nuclei: 

Tr^i = (TM]T< M q M -q>- = (T/^)]>>2 q 2 D(q, W ). (8) 
q q 

Here a x = a± or ay, depending on the direction of the magnetic field. The a's are undetermined 

numerical factors which depend, besides the couplings in Eq. (3) on details of the structure and 

chemistry. In (8) q is two-dimensional; the magnetization correlations between the planes is purely 

from the very small long-range electro-magnetic effects, and assumed smaller than the temperature. 

Now we note that near x = x c (ignoring log T corrections) 
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Um(u -► 0)(T/w) q 2 ImD(q, u] ' 



70 (T) 



+ 0(1/E F ) 



(9) 



Since 70 (T) = r ~ 1 + AT at G=0, a scale-invariant contribution to the local fluctuation on Cu indeed 
arises in the pure limit, AT » t^ 1 . 

We write the ordinary contribution to nuclear relaxation rate of a nucleus n with the field 
applied in the direction i as 



^ = lim^oEr F i( ( li)lVo / (q^)/Ml- (10) 



nr-p — 1 

q,i 

In (10), n Fi(q) include the measured hyperfine constants as well as the lattice form factors, as 
given, for instance, in Ref. (18). i are orthogonal to i. To this we supplement for Cu the anomalous 
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contribution coming from (9), 

Cu T -i = T 

lan ' ; W s7 o(T)' 1 ) 

Equation (10) has been evaluated extensively, with x'o calculated in various approximations 19 ' 12 
and with the form of x'oili u ) deduced from experiments 9,10 for Lai.ggSr 14CUO4 which show in- 
commensurate antiferromagnetic correlations. 10 We defer the role of such correlations, which in any 
case are of an ignorable magnitude in 12306.93, an d write (10) in the Fermidiquid form, so that 

Cllrp-1 _ XO rnlCu A_|2 1 a i 1 /io\ 

Tl ' i_ ^ 1 11 + wI^(T)' (12) 

and 

° T ^(i) T W < 13 > 

where in n Ai, we have absorbed the hyperfine constants, form factors and other factors after scaling 
out a factor (xo/Vb-Ef)- Note also that oi\ has the same dimensions as A{. 

We now quantitatively compare (12) with the experimental results. We first note that lim q — > 0, 
to — ► Re q 2 D(q, oj) = a 2 /W s , so that for a 2 <C 1, this provides a negligible contribution to the 
uniform magnetic susceptibility and the Knight shift. So the hyperfine constants deduced from the 
measured Knight shift should be used only in the usual (Fermidiquiddike) contributions to T^ 1 . 
At high temperatures, T » r^ 1 , (12) behaves as a + bT; the constant part is replaced by the term 
proportional to T tq for T » t^ 1 . If three dimensional coupling of the anomalous magnetization 
fluctuations were significant, it would also produce a term proportional to T. 

In Fig. (2) we compare the calculated result for Cu T]~^ L from Eq. (12) with field in the plane with 
the experimental results in (12306. 93). We deduce | Cu A^| 2 by scaling up the measured (TiT) -1 
of oxygen by the ratio of the appropriate measured hyperfine coupling constants. We can estimate 
70 (T) from the measured resistivity. The fit in Fig. (2) is with A « 1 and Tq 1 w 25K, similar 
to what one deduces from the measured resistivity of (12306.93 )• This leaves the undetermined 
parameter, a]_/W s , to fit the data. The fit in Fig. (2) is with «j_/W s 3.2(msec)~ 1 , which with 
W s ~ Ep ~ 10 4 K implies |a±/ Cu A ± / 2 w 6 x 10~ 2 . 



We also compare the Cu Ti measured by NQR in a single crystal of Lai.seSr 14CUO4. The 
only significant difference in parameters needed to fit the data is a larger value Tq" 1 10 2 K. This 
is consistent with the fact that the extrapolated residual resistivity of optimum (124) crystals is 
generally higher than that of optimum (123) crystals. The fit in Fig. (2) uses aj_/W s « 2.8(msec)~ 1 . 

To see if a\ is reasonable, we compare it to a hypothetical situation where d x 2_ y 2 and d xy 
orbitals on Cu are degenerate, so that a local orbital as well as a dipolar fluctuation contribution to 
T^f 1 exists with field in the plane. The ratio |A or b/Aj_| 2 in this case (where both relaxation rates 
are ~ T) can be easily calculated 20 and is ~ 20. We have two reduction factors in a 2 compared to 
A 2 rb , that due to ((j) 2 } ~ x, and more importantly, because the distance to the neighboring oxygen 
p-orbitals is involved. A reduction factor of O(10 -2 ) is therefore not unreasonable. 

In general (aj_/a||) 2 > 1 is expected, the details depending on structure. Experimentally this 
ratio is about 2.6 in Lai.ssSr 15CUO4 and about 4 in YBa2Cu306.g. This is in the right direction 
from the consideration here since the former is a single layer and the latter a double layer compound. 
Measurements of anisotropy in other compounds appear not to have been made. We would predict 
that this ratio rises in compounds with larger number of layers. Any numerical estimates do not 
appear feasible. 

At compositions away from the ideal, i.e. go|x — x c | >> T/To, the fluctuation modes of eq. (1) 
have a gap. This has the same effect in T^ 1 as an enhanced t^ 1 . This is qualitatively consistent 
with the experimental results that there is a smoother variation of T^f 1 from low temperatures 
to high away from the ideal composition. It is worth noting also that a finite G also leads to a 
cross-over at temperatures of 0(G) from a lower effective mass to a higher effective mass. This 
should be observable in thermodynamic properties: both specific heat and magnetic susceptibility 
as indeed it is. 21 

As has been discussed before, 12 ' 19 the °T^ 1 in (I23O7) is consistent with (13) and given quan- 
titatively in terms of the measured xo an d hyperfine constants. This is true in the present theory 
at all x. For °T 1 _1 in Lai.ssSr 15CUO4 to be consistent with the neutron scattering results, the 
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absolute magnitude of the q-dependent part of x"(q, oj) at low to should be about a factor of 3 less 
than used in Ref. (10) or the deduced (hyperfine constants) 2 smaller by a similar factor. 

In contrast to T^f 1 which is related to fluctuations integrated over q, T^ 1 of Cu, as derived 
by Pennington and Slichter 2 is dependent primarily to fluctuation at q « 7r/Rc u -Cu- So anti- 
ferromagnetic correlations, however weak are picked up in T^ 1 . Berthier et al. 11 find using the 
neutron scattering results as a function of temperature in 12306.7 and in 12306.93 that the mea- 
sured correlation length of only about two lattice constants in the former and one or less in the 
latter accounts for the measured T^" 1 (T). Similar conclusions have been arrived at by Walstedt for 
Lai.86Sr 14CUO4. 

^From Eqs. (1), (4) and (7) one can deduce that the orbital magnetic fluctuations contribute 

ImXorb(q,w) - Hb (^j (a 2 q 2 )ImD(q,cj) (14) 

to the total magnetic fluctuations, Imx(q, lj), with a form factor such that they are seen only at 
the Cu (or symmetry equivalent) sites and absent at the oxygen (or symmetry equivalent) sites. 
Here a ~ 1.9A is the nearest neighbor Cu-0 distance, and ao is of the order of the atomic radius. 
To this should be added the usual spin-fluctuation contributions. Should there be significant (but 
not singular) AFM correlations, for example due to nesting, the Fermi-liquid contributions to Eqs. 
(12) and (13) can easily be modified to include their effect. 

A contribution, smooth in both q and w extending up to high energies, as in Eq. (14) is indeed 
observed. 22 To test Eq. (14) in detail requires measurements up to high frequencies at various 
temperatures. The only measurement reported 22 up to high frequencies has been done only at 
T = 17K in Lai.gsSr 15CUO4. Fig. (3) of Ref. (22) gives the frequency dependence of the integrated 
absorption as very slowly decreasing to ~ 0.2 eV below a weak peak at ~ 20 meV. This is quite 
consistent with Eq. (14). Measurements at other temperatures in that compound and up to high 
frequencies and various temperatures in 12306.9 are urged. Especially crucial are measurements of 
the fluctuation spectrum in several Brillouin zones to deduce separately the fluctuations centered 
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on Cu alone and on O alone The former should see (14) while the latter should be the usual Fermi- 
liquid form, i.e. « N(0)w/qVF for uj < qVp and f» beyond. These are hard experiments, but they 
should resolve the principal problem in the field. 

I am pleased to acknowledge very helpful discussions with A. Sengupta, G. Kotliar, P. B. 
Littlewood and Q. Si about this work, and extensive discussions on NMR with R. E. Walstedt, H. 
Alloul, C. Berthier and C. Hammel and on neutron scattering with G. Aeppli. 



10 



References 

1. For a review, see B. Batlogg in High Temperature Superconductivity, ed. by K. Bedell et al., 
Addison Wesley, New York (1989). 

2. For review, see C. H. Pennington and C. P. Slichter in Physical Properties of High Temperature 
Superconductors II, D. M. Ginsburg, editor, World Scientific, Singapore (1989). 

3. P. C. Hammel et al, Phys. Rev. Lett. 63, 1992 (1989) 

4. H. Alloul, T. Ohno and P. Mendels, Phys. Rev. Lett. 63, 1700 (1989) 

5. R. E. Walstedt et al., Phy. Rev. B 38, 9299 (1988). 

6. C. M. Varma et al., Phys. Rev. Lett. 63, 1996 (1989). 

7. B. S. Shastry, Phys. Rev. Lett. 63, 1288 (1989). 

8. J. Rossat-Mignod et al., Physica 192 B, 109 (1993); J. M. Tranquada et al. Phys. Rev. B4.6, 
5561 (1992); H. A. Mook et al., Phys. Rev. Lett. 70, 3490 (1993); L. P. Regnault et al., 
Physica C 235, 59 (1994); B. Keimer et al., (preprint). 

9. T. E. Mason et al., Phys. Rev. Lett. 71, 919 (1993). 

10. R. E. Walstedt, B. S. Shastry and S.-W. Cheong, Phys. Rev. Lett. 72, 3610 (1994). It is 
also worth noting that the result reported (T. Imai, Phys. Rev. Lett. 70, 1002 (1993)) that 
at high temperatures Lai.ssSr 15CUO4 has the same Cu T 1 _1 as La2SrCuC>4 has been modified 
by more recent measurements (M. Matsumra et al., J. Phys. Soc. Japan 65, 699 (1996)). 

11. C. Berthier et al., Physica C 185-9, 1141 (1991); Physica Scripta T49, 131 (1993) and private 
communication; R. Walstedt (unpublished). 

12. P. B. Littlewood et al., Phys. Rev. B 4.8, 487 (1993). 

11 



13. C. M. Varma, S. Schmitt-Rink and E. Abrahams, Solid State Comm. 62, 681 (1987). 

14. C. M. Varma, Phys. Rev. Letters. 75, 898 (1995) and Preprint (1996), submitted to Physical 
Review and available on http: / /www.lanl.govj as |cond-mat / 96071 05 , In this preprint, some 



important modifications of the results in the PRL are given. In particular, a damping mecha- 
nism of the collective mode due to particle-hole scattering near the Fermi-surface leads to the 
form of T) in Eq. (1) which is more important than Landau damping for all momenta. 

15. In the tight-binding approximation the phase factors for electronic wave- functions is given 
only at the lattice points. The results derived below depend on varying the phase of the 
wave-functions between lattice points. Therefore Bloch waves must be used for the electronic 
wave- functions. 

16. P. A. Lee and N. Nagaosa, Phys. Rev. 43, 1223 (1991). 

17. For general k's there is also a contribution from the first term in (5) with properties similar 
to those evaluated here. Since we are unable to evaluate any of the matrix elements quan- 
titatively and the aim is only to show the feasibility of the idea and fit to experiment with 
reasonable parameters, we do not discuss these. 

18. (a) F. Mila and T. M. Rice, Phys. Rev. B 40, 11382 (1989); (b) A. J. Millis and H. Monien, 
Phys. Rev. B 45, 3059 (1992). 

19. Q. Si, Y. Zha, K. Levin and J. P. Lu, Phys. Rev. ^7, 9055 (1993) and references therein. 

20. Y. Obata, J. Phys. Soc. Japan 18, 1020 (1963). 

21. J. W. Loram et al, Physica C 235-240, 134 (1994). 

22. S. M. Hayden et al., Phys. Rev. Lett., 76, 1344 (1996). 



12 



Figure Captions 



Figure 1: Illustrates the physical content of Eq. (10). A Cu d x 2_ y 2 orbital is surrounded by the 
p x ,y orbitals of oxygen atoms with the phase as shown. A deviation of the phase ~ q in the direction 
shown produces a fluctuation with d xy symmetry (shown shaded) about the Cu-site leading to a 
orbital magnetic moment fluctuation ~ q. 

Figure 2: Experimental results from Ref. (5) for Cu nuclear relaxation rate with field the 
planes in YBa2Cus06.9, and nuclear quadrupole relaxation rate in Lai.ssSr 15CUO4 from Ref. (10) 
compared with the calculations in this paper. The two parameters required are discussed in the 
text. 
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